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INTRODUCTION 


On the national forests there has been a trend during recent years 
toward placing on a systematic basis such fire- control practices as 
rating fire danger, determining the proper size of suppression crews 
and speed of attack, and planning fire-suppression strategy. This 
trend has resulted in greatly improving the techniques of planning and 
managing the fire-control organization. At the same time it has 
revealed a serious lack of essential information. The outstanding 
need is for specific data on the rate at which forest fires may be 
expected to spread under various conditions of forest cover, weather, 
and topography. 

Earlier studies in California by Show (1/3) * and Curry and Fons 
(6, 7) consisted of observations on the rate of fire spread under natural 
forest conditions of fuel, weather, and topography. The initial series 
of experiments by Curry and Fons (6) in one of the least complex fuel 
types, ponderosa pine needles, provided data on the empirical relation- 
ship between the rate of fire-perimeter increase, litter moisture content, 
wind velocity, slope, and time from origin of fire. Similar observations 
were made on other fuel types. 

Analysis of these data made it apparent that without a basic under- 
standing of the many physical processes involved, the results of the 
experiments could not be applied beyond the limits of the conditions 
under which the observations were made. Therefore, it became neces- 
sary to determine the fundamental laws governing the rate of fire 
spread in forest-type fuels. Under field conditions, however, none of 
the important factors, such as the attributes of the atmosphere, the 
arrangement of the fuel bed, and the physical properties of the fuel 
particles, remained sufficiently uniform throughout an experiment to 
allow exact description of the numerous variables influencing the rate 
of fire spread. In order to understand the influence of these variables 
it was decided to conduct experiments with model fires in beds of 
chosen homogeneous fuel particles. These were at first performed by 
Curry and Fons (7) in still air for different conditions and kinds of 
fuel, fuel-moisture content, and fuel-bed compactness, and later by 
Fons (9) in a wind tunnel in which air velocity was controlled. 
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Concurrently with the latter experiments, an analytical treatment 
of the rate of spread of fire was being developed and utilized to direct 
the course and improve the techniques of fire control. This experiment 


and its results are described in this paper.® 
GLOSSARY OF SYMBOLS 


A=cross-sectional area of fuel particle, sq. in. 
A,=surface area of particle, sq. in. 


a, wg = . . .. n=subscripts designating different fuels. 
W, W. Wa 
wee to 
=A: Avys iat ; AnYn_ 


C=proportionality constant relating a particular kind of fuel bed to a standard 
chosen bed. 


C,=specific heat of dry fuel at constant pressure, B. t. u./lb. ° F. 
C,=specific heat of moist fuel at constant pressure, B. t. u./Ib. ° F 
C,=temperature-difference ratio coefficient. 
D=diameter of fuel particle, inches. 
d=prefix indicating differential, also depth of fuel bed, inches 
E=fuel particle shape factor, oA”. 
f=over-all transfer factor, also suhseript indicating flame or film. 
f-=film conductance for convection, and conduction, B. t. u./sq. ft. hr. ° F 
f,=heat transfer factor for radiation, B. t. u./sq ft. hr. ° F. 
h=heat of adsorption of dry fuel, B. t. u./lb. 
ky=thermal conductivity of the air film, B. t. u./sq. ft. hr. ° F./ft. 
.=spacing between particles, inches. 


M=moisture content, weight of water/weight of dry fuel. 
N=number of particles. 


q=total rate of heat transferred, B. t. u./hr. 


qe=rate of heat transferred by convention and conduction, B. t. u./hr 
qr=rate of heat transferred by radiation, B. t. u./hr. 


R=rate of spread or e in r-direction, ft./hr. 
S=total fuel surface in volume Vj, sq. in 
T=temperature, degrees Rankine. 
t=temperature of fuel particle at time 6, ° F. 
to=temperature of fuel particle far away from flame, ° F 
t;=temperature of fuel particle near flame, ° F. 
i;=temperature of flame, °F. 
i;=ignition temperature of fuel particle, ° F. 
V=wind velocity at 1 foot above fuel-bed surface, m. p. h. 
V,=volume of fuel bed, cu. in. 
V.=volume of fuel in fuel bed volume Vj, cu. in 
v=volume of particle, cu. in. 
W =weight of dry fuel, lbs. 
= (gamma) density of dry fuel, lbs./cu. ft. 
y= (gamma bar) density of moist fuel, lbs./cu. ft. 
e= (epsilon) flame emissivity. 
@= (theta) time in hours. 
6;=ignition time of fuel; time in hours to raise a fuel particle from its initial 
temperature to ignition temperature. 
= (lambda) volume of voids per unit of fuel surface, inches. 
v¢=(nu) kinematic viscosity of the film, sq. ft./hr. 


o=(sigma) surface-volume ratio, or fineness of fuel, inches~! (1--number of 
inches). 


= (phi) slope of terrain, percent. 


5 A portion of the analysis in its earlier stages was submitted by the author in 
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DERIVATION OF EQUATIONS FOR RATE OF SPREAD IN AN 
IDEALIZED HOMOGENEOUS FUEL BED 


Forest cover may be regarded as a heterogeneous mixture of dead 
and green fuel particles of varying size, shape, and density. Two 
layers may be distinguished: a lower layer of inflammable litter next 
to the ground surface, and an upper layer of dead and green standing 
vegetation. A fire may spread through either or both layers at a 
definite rate depending on the magnitude of many variables pertaining 
to (1) the atmosphere, (2) the fuel bed, (3) the fuel particles, and (4) 
the topography. The fire will spread quickest in that direction and 
through that layer where ignition is most favorable, or where resist- 
ance to flame propagation is at a minimum. 

A forest fire characteristically travels more or less rapidly through 
the inflammable cover, singeing superficially only the finer, more 
ignitible particles, later slowly and thoroughly consuming the great 
mass of heavier fuels. In this theory of the rate of spread it is im- 
portant to center our attention on the action at the head of the fire, 
where there is generally ample oxygen to support combustion. At 








PARTICLE NUMBER 2-2 T-i Ft Nei 
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Figure 1.—Graphic representation of the rate of spread of fire in an idealized 
fuel bed composed of identical rodlike fuel particles. 





7-0 P- DIRECTION <> v= @ 


this point only sufficient heat need be produced for a self-supporting 
fire to bring the immediate adjoining fuel to ignition temperature. 
Accordingly, flame propagation, or rate of fire spread in a fuel bed, 
may be visualized as proceeding by a series of successive ignitions, its 
magnitude controlled primarily by the ignition time of the particles 
and the distance between them. 

As in other complex natural phenomena, the problem must be 
idealized in order to facilitate mathematical treatment of the theory: 
It is believed, however, that doubt as to the applicability to field 
conditions of the theoretical equations developed by this idealization 
will be removed when the predicted effects of the individual variables 
on rate of spread are compared with previous experimental results 
and with observations under actual field conditions. 


RATE OF SPREAD AS A FUNCTION OF PARTICLE SPACING AND IGNITION TIME 


For this derivation an idealized bed is assumed (fig. 1) composed 
of identical rodlike fuel particles, standing upright at equal distances 
from each other. If 6; is the time for the n particle to burst into 
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flame after the (n—-1) particle has ignited, and if Z is the distance 
between the particles, it follows that the rate of spread R is 


R= (1) 


dr . ie j 
where R= 4 in the r-direction down the row of fuel particles. 


RATE OF SPREAD AS A FUNCTION OF FUNDAMENTAL VARIABLES 


The ignition time 6; in equation (1) depends on two types of varia- 
bles: Those dealing with (a) the fuel particle and (b) the atmosphere 
surrounding the particle. 

Again it is assumed that the bed is composed of idealized fuel 
particles (fig. 1) spaced at a uniform distance Z and that immediately 
upon ignition of the (n~1) particle its flame at temperature t, extends 
far enough to be in contact with the n particle which is at a tempera- 
ture ¢t,. When the latter particle is exposed to the gases of flame 
temperature, t,, it will absorb heat by convection, conduction, and 
radiation. 

In the theory of fluid flow (11, 14) it is generally accepted that when 
a gas or liquid moves over a solid a thin film of fluid forms on the 
surface. This becomes thinner as the velocity of the fluid increases, 
and the thinner the film the more rapid will be the transfer of heat to 
the solid by conduction, convection, and radiation. The film con- 
ductance or transfer factor, f, is used to express the net rate of heat 
transfer through the film per unit area per degree difference in tempera- 
ture. Factor f is the sum of two transfer factors, f, for conduction 
and convection, and f, for radiation. These are determined inde- 
pendently because the heat exchange by radiation depends primarily 
on the temperatures of the surfaces, while that for conduction and 
convection depends on variables such as the physical properties of the 
fluid, the shape of the solid surface, and the velocity of the fluid past 
the solid boundary. The method of calculating the transfer factors 
is discussed below. 

The heat transferred by convection and conduction in unit time 
from the atmosphere at t,; to the particle is 


Qe=f-As(ts—1) (2) 


The heat transferred by radiation in unit time from the atmosphere 
at t, to the particle is approximately: 


r=, As(ts— t) (3) 
The total heat absorbed in unit time is 
q=qet+¢qr=(fe+f,) (t;—t)A,g (4) 


where 
A,=surface area of the particle; 
t;= uniform temperature of flame surrounding the particle; 
t=temperature of particle at any time 6, assumed to be uniform 
throughout the volume of the particle, since the temperatuce 
gradient is very low for fine fuel particles with large surfare- 
volume ratio; 
f, and f, are the coefficients of heat transfer for convection and radia- 
tion energy, respectively. 
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The rate of heat absorption, qg, will cause a rate of temperature 
dt , : a 
change 7’ according to the density of the moist particle, 7, the 


specific heat of the moist particle, C,, and the volume of the fuel 
particle, v, as follows: 


dt ‘ 
q=C, 6 (5) 
Substituting (5) in (4) yields: 
_ dt _ (fe- +f )As (6) 
(t;—t) yyw 
By integrating equation (6), with Geta * assumed constant, 
wf 10 


equation (7) is obtained: 


~in(t, 9 =e 
yCw 


where /n represents the Napierian logarithm. 


+lnC (7) 


The initial conditions to be satisfied are: when 6=0, t=t,. 
By substituting these values in (7) the constant of integration becomes 


lnC= —In (t;—t,) (8) 
Substituting (8) in (7) and rearranging 
_ Ce +- f,) 08 
* ee 7C, ) 


c=A, : 
where — ° “=the surface-volume ratio of fuel particle. 


The final conditions to be satisfied are: when 6=6,, t=t,, where 
6,=ignition time, and t;=ignition temperature. 

Making these substitutions and rearranging for 6,;, equation (9) 
becomes 


0,= VC, In =a) (10) 


o(fettr) Be 


Equation (10) is an expression for the ignition time, 6;, in terms 
of the physical properties of the light fuel particle,® the flame tempera- 
ture, and the heat transfer factors. 

The rate of spread for an idealized homogeneous fuel bed is derived 
directly by substituting the equivalent of 6; from equation (10) in 


equation (1): 3 
RaaLethe (11) 


v0, i) 
In Ss 7 


In this equation R is expressed in terms of certain fundamental 
variables that control the spread of a fire in light forest fuels. Such 
6 Derivation of an expression for heavier fuel particles, in which there is an 


appreciable temperature gradient through the particle, is possible by the applica- 
tion of an appropriate form of the conduction equations (4). 
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quantities as wind velocity, moisture content, time, and topography, 
which are ordinarily regarded as controlling spread, do not appear in 
the equation. They affect the rate of spread indirectly through 
their influence on the more fundamental variables. 


Firm ConDUCTANCE FOR CONVECTION 


Experimental data for heat transfer from gases to single cylindrical 
rods, with the gas flowing at right angles to the axis, have been 
collected by numerous investigators. Data for cylinders from 0.0004 
to 3.75 inches in diameter (1/1, p. 218), applicable to fine forest fuels 
such as needles and twigs, have been summarized and correlated in 
terms of dimensionless quantities related by the following empirical 
equation, which may be used as an approximation for calculating f.: 


7 \ 0.56 
LP -0.45+0.33(*") (12) 


where f, =film conductance for forced convection 
D=diameter of cylinder 
k,=thermal conductivity of the air film 
v-=kinematic viscosity of the air film 
V=air velocity. 





Thermal conductivity and kinematic viscosity depend on the 
arithmetic mean temperature of the film. Values of these properties 
for air are given in the International Critical Tables for a wide range 
of temperatures. 

Forest fuels such as needles and twigs are nearly enough cylindrical 


: ‘ ; ; ae , 
in shape to allow replacing D in equation (12) with i where o is the 


surface-volume ratio of the fuel particle. 
Hpat TRANSFER COEFFICIENT FOR RADIATION AND FLAME TEMPERATURE 


The net rate of heat transfer by radiation, g,, to a small particle 
surrounded by a flame (14) is obtained by the following equation: 


Seats 7, eo Le 4 
=0.172¢A,| (54) Gas) | (13) 


where e=flame emissivity ’ 
T,=flame temperature 
1=particle temperature 
A,=surface area of the particle 
An expression for f, may be derived by substituting simplified 
equation (3) in general equation (13), 


; 6) - i) 
‘ 0.172¢{ 100 (06) | 


a ee am 





7 Radiant emissivity is defined as the ratio of the total emissive power (total 
radiant energy emitted per unit time from unit surface) of a surface to the total 
emissive power of a black body at the same temperature. 
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Flame temperatures were measured by chromel-alumel thermo- 
couples held at several points in the luminous flames of different kinds 
of burning forest fuel. The average temperature so measured was 
approximately 1500° F. Simultaneously, the emissive power of the 
flames was measured by a portable radiometer described by Boelter 
(2). From the measured flame temperature, the radiometer element 
temperature, and the flame emissive power, ¢,, the emissivity, ¢, was 
calculated by the radiation equation (13). The approximate value 
of « so determined was 0.227. 

Until more data are acquired, the values used for ¢ and ¢; in calcu- 
lating f, are 0.23 and 1500° F., respectively. Since the particle is 
being heated until it reaches ignition temperature, an approximation 
must be made for the fuel temperature, ¢;. On the assumption that 
the rate of change of the particle temperature is constant, the fuel 
temperature, t,, may be approximated by using the arithmetic mean 
of the ignition temperature, ¢;, and the fuel temperature, ft). While 
t, may vary over a considerable range between ¢, and fo, its approxima- 
tion at the mean value of these two will not result in an appreciable 
error in the calculation of f,. 


IGNITION TEMPERATURE 


There is no general agreement as to the ignition temperatures of 
solid fuels. The values determined by several investigators whose 
work was reviewed by Brown (3) show a large variation. For 
instance, values as high as 1200° and as low as 650° F. have been 
obtained for sawdust and as low as 410° for wood shavings. The prin- 
cipal causes of the variations are (1) the definition of ignition temper- 
ature, which may refer to the specimen temperature or to the external 
temperature of a bar, an ambient gas, or an oven, and (2) the criteria 
chosen by different investigators to indicate the ignition point. 
Criteria that have been used are the appearance of glow, the appear- 
ance of flame, or some critical point in a temperature- or pressure- 
time curve. 

To be consistent with the preceding theory and the equations 
developed for the rate of fire spread in light forest fuels, the following 
definition is proposed: The ignition temperature of a light forest-fuel 
particle surrounded by ambient hot gases is the average specimen 
temperature at the instant a flame appears. For light forest fuels of 
known physical properties exposed to hot gases of known tempera- 
ture, the average specimen temperature may be calculated by equa- 
tion (10) from the ignition time, which is determined experimentally. 
Until such experimental data are available, an arbitrary chosen value 
of 550° F. is recommended for all light forest fuels irrespective of fuel 
size and moisture content. Jones and Scott® found the ignition 
temperature of dry wood to be in the range of 518° to 554° F. 


ParRTICLE SPACING AND SURFACE-VOLUME RaTIo 


To calculate the rate of fire spread in a fuel bed by equation (11), 
L, the spacing, must be related to the measurable fuel and fuel-bed 
variables. 

8 Jones, G. W., and Scott, G. S. CHEMICAL CONSIDERATIONS RELATING TO 


FIRES IN ANTHRACITE REFUSE. U. 8S. Bur. Mines Rpt. Invest. 3468, 13 pp., 
illus. 1939. [Processed.] 
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For a bed of fuel, the ratio of voids to fuel-surface area is 


i Seat: 


= yo 15) 
S ( 
where \=volume of voids per fuel-surface area 
V,=volume of fuel bed 
V,=volume occupied by the fuel in volume V, 
Y . y 
S=total surface area of fuel in volume V,. 
The reciprocal of \ may be thought of as the compactness of the fuel. 
The volume of the idealized fuel bed is 


V,=L’dN (16) 


where L=spacing of particles 
d=depth of fuel bed 
N= total number of particles. 
The volume occupied by the fuel in the bed of volume V, is 
, W a 
V.= (17) 
sf 


where W is the total weight and y is the density of the dry fuel. The 
total surface area of the fuel in bed of volume V;, is 


oW 
é | 


S = 


Substituting equations (16), (17), and (18) in (15) gives 


L?dN—W/+y 


a (19) 


Solving equation (19) for Z and letting the cross-sectional area 
om \. ‘ ; : ’ 
A= Nvd (since for an idealized bed the length of the particles is equal 
to the depth of the bed) (19) reduces to: 


L=A"(1+on)" (20) 


Equation (20) expresses the particle spacing in terms of the cross- 
sectional area, the surface-volume ratio of the fuel particle, and the 
ratio of voids to fuel-surface area of the bed. 
Multiplying both sides of equation (20) by o and letting H=ocA” 
then 

oL= E(1+ oh)” (21) 


The quantity F is independent of the size of the fuel particle and 
depends only on the shape of the cross section. For instance, for 
round sticks regardless of size, HE=22; for square sticks, H=4.0; 
for more irregular cross sections, E>4.0. The cross-sectional area 
of an irregularly shaped particle may be determined from the relation 
A=, where p is the perimeter and o the surface-volume ratio. 
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Speciric Heat 


Since the fuel particles nearly always carry adsorbed water, the 
specific heat of the moist fuel, C,, may be expressed in terms of 
specific heat of dry fuel, its moisture content, and the heat of adsorp- 
tion. 

The total heat required to raise a moist fuel particle at temperature 
tp to ignition temperature ft, is 


H=(W+MW)C,(t,;—t) (22) 


where W=dry weight of fuel 
M=moisture content, weight of water per unit weight of dry 
ms fuel 
C,=specific heat of moist fuel. 

The total heat, H, has three components: the heat required (1) to 
separate the water from the fuel, that is, the heat of adsorption; (2) 
to raise the dry fuel from temperature ft) to ignition temperature, t,; 
and (3) to vaporize the water and superheat the vapor to temperature 


“The heat required to raise the dry fuel to ignition temperature 1s 
A,= WC, (ti—to) 23) 
where C, is specific heat of dry fuels (8)°; 
and to raise the water to temperature ¢; is 
H,= MW (hi—ho) 


where h;=1,057 +0.46t,, the heat (/, p. 309) of the superheated steam 
at low pressure in B. t. u. per pound, and hp=t)—32, the heat content 


of water at temperature t,, whence 
H.=MW (1,089 +-0.46t,—ty) (24) 


The heat of adsorption per unit dry weight of fuel is h, and heat 
required for weight of dry fuel W is 


H,=Wh (25) 


The heat of adsorption for wood at 32° F. is approximately 32 B. t. u. 
per pound of dry fuel (10, p. 286). 

Adding equations (23), (24), and (25), the total heat, 7, to bring 
W, weight of moist fuel, to ignition temperature ¢; is 


H=WC,(ti—to) + MW (1,089 +-0.46t;— ty) + Wh (26) 


®The equation for specific heat of dry wood from results considered as prelimi- 
nary is 

: i YU fis ‘) 

Cp=0.246 + 0.00063 ( 2 B.t.u./Ib. °F, 


679977—46 2 
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Equating (26) to (22) and solving for C,, the specific heat of moist 


fuel is 
= 1 ,089 +0. 46t;— — to 
Coat gg | (Oe ot t+, | (27) 


Density or Moist Fue. 


The relation between the density of the moist fuel and dry fuel, 
y, is y=(1+M)y. Ordinarily the density of the dry fuel, y, and 


CF.) 
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DIRECTION OF FIRE SPREAD > 


Figure 2.—Diagrammatic fuel temperature gradients at and near the flame front: 
B is for a fire burning with a higher wind velocity or on a steeper slope than A. 





the moisture content, M, are obtained experimentally, and ¥ for the 
moist fuel is calculated. 


Fue TEMPERATURE 
Temperature ¢, (fig. 2) will obviously vary with the following: 
(1) The flame angle from the horizontal, which is controlled by the 
wind velocity and slope of terrain; (2) the size of the flame front, 
which is controlled primarily by the time from the start or “age of 
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fire’ and the amount of fuel burning at the head of the fire. Tem- 
perature ¢, varies with the flame angle because, as the angle becomes 
smaller, more hot gases pass directly over the unburned fuel bed in 
the direction of the main flame front; thus the fuel particles take on 
more heat by convection and radiation. Temperature ¢, varies with 
the size of flame, which governs the amount of heat transferred to the 
unburned fuel particles by radiation. The time from start of fire also 
affects temperature ¢,, because at zero time ¢, must equal ft, but at 
subsequent intervals the fuel particles ahead of the fire have taken 
on heat, causing ¢, to become greater than 4. Temperature ¢, is not 
easily measured even under ideal conditions. Temperature t) of fuel 
far removed from the advancing flame front is, however, readily 
determined. Temperature ¢, may be eliminated by substituting ty 
for it and correcting this new temperature-difference ratio in equation 
(11) by a coefficient, C,, in the following manner: 


t,—t t,—t 
Let rik pe °) stab | 28 
: ar bk (28) 
wherefore C, == ; (29) 
0 


Figure 2 shows diagrammatically the fuel-temperature gradient 
and temperature differences’ at and near the flame front. For all 
practical forest conditions ty; may be assumed to be about 1500° F. 
The effect of increasing wind or slope makes t, numerically larger, 
which makes C, smaller, according to equation (29). Since ¢, in- 
creases numerically with time from start of fire and size of flame, 
these also will make C, smaller. Any independent change in f, will 
be nearly compensated by a change in the same direction in ¢, and will 
not therefore greatly change C,. The temperature-difference ratio 
coefficient is primarily a function of wind velocity, slope, size of flame, 
and time from start of fire and must be evaluated experimentally. 
The approximate numerical range of C, may be determined by substi- 
tuting a reasonable value of ¢ and the limiting values of t, in equation 
(29). For the field conditions of tj=80° F. and t;=550° F., C; will 
vary from 0.67 to more than 1.00. At the very start of a fire, when 


j tj=t, C,=1.00. For the extreme conditions of zero flame angle, 
‘ which occur at extremely high wind velocity or infinite slope, C, will 
approach 0.67 when ¢, approaches ¢;. For fires burning in still air in 


dead fuel exposed to full sunlight, C, may become somewhat greater 
than 1.00 with time from start of fire; for under these conditions t, 
will be high and the fuel may actually be cooled by the indraft. 


TRANSFORMATION OF EQUATIONS DERIVED FOR IDEALIZED 
HOMOGENEOUS BEDS TO NATURAL HOMOGENEOUS BEDS 


In the idealized fuel bed all particles were assumed to be vertically 
oriented with equal spacing. In the natural homogenous bed, such 
as a mass of fallen pine needles, more or less horizontally oriented at 
random, the spacing is not equal but has instead voids of varying 
sizes. There are other differences between idealized and natural 
beds. In the latter the particles are in occasional contact with one 
another, and will not always be oriented with the flow of hot gases as 
in the idealized bed; this will influence the rate of heat conducted into 
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7 0 There may also be a fuel-moisture gradient in the natural 
uel bed. 

Granting that all the variables affecting the rate of spread have 
been included in the derivation, according to the theory of similitude 
(12, p. 78), if the processes of burning in the two beds follow the same 
laws and the respective variables stand in a constant ratio to each 
other over their whole range, the application of equation (11), derived 
for idealized beds, to fires in natural beds is permissible. On the 
assumption that the individual variables in the two processes are 
respectively proportional, as required by the theory of similitude, 
equations relating these variables may be written; for example, the 
relation between spacing L, of the idealized and L, of the natural bed 
may be expressed as 


L, — CL, 


where C; is a constant quantity for the whole range of spacing. Other 
variables such as specific heat, density, etc., may be dealt with in a 
similar manner. 

Accordingiy, the rate-of-spread equation will be identical in the 
natural and idealized beds, except that each variable will be multiplied 


‘ : ss = — —f 
by its corresponding coefficient. With replaced by oT ), 
seme | | ee | 


in which C, serves as both the temperature-difference ratio coefficient 3 
and as a coefficient relating the natural to the idealized values of 
ty and t, equation (11) for a natural fuel bed becomes 


pa CMfet CrCl 


C,5CC,0 snl (E *)) (30) 


Inasmuch as the coefficients C;, C,, etc., remain constant 





let /= i me (31) 
CEAS) = 


oi a snl C, = — ) es 


where C represents a proportionality constant relating the rate in a { 
natural bed with the rate in a standard bed under identical conditions. 

The coefficients C, and C in equation (32) are evaluated experi- 
mentally. As shown in the above discussion, C; is primarily a func- 
tion of wind velocity, slope, size of flame, and time from start of fire. 
The numerical value of C is dependent in part on the type of fuel 
bed, that is, whether it is a mass of fallen litter or standing vegeta- 
tion, such as grass or brush, and in part on the technique used in 
measuring the variables, such as wind, fuel moisture, and fuel tem- 
perature. To determine the numerical value of C, for any particular 
wind, time, and slope it is first necessary to evaluate C for the par- 
ticular fuel used in the experiments. However, in order to obtain 
the value of C with reasonable accuracy we must have data for a wide 
range of fuel temperature, and this is experimentally impracticable. 
On the other hand, if C, be evaluated with C taken as unity for a 


then 
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chosen standard natural fuel bed, for which a precise technique for 
measuring the wind, moisture content, and other variables has been 
developed (to a point where measurements of all variables can be 
duplicated), then C can be relatively evaluated for any other kind of 
natural bed. This is done by determining R experimentally for a 
particular kind of bed and computing a value for C by equation (32), 
using C, as evaluated for the standard bed. The coefficient C deter- 
mined in this manner would then be a constant, relating the particu- 
lar kind of natural bed not with the idealized bed, but with the chosen 
natural standard. 


DERIVATION OF EQUATIONS FOR RATE OF SPREAD IN NATURAL 
HETEROGENEOUS FUEL BEDS 


Fuel beds of the homogeneous type are of rare occurrence in nature. 
Normally a natural bed is composed of particles varying in size, shape, 
density, and specific heat. Since the ignition time, ¢; is dependent 
on the material and dimensions of the particles, it follows that the 
rate of spread will be different for different particles. 

Consider an idealized heterogeneous bed composed of particle 


types a, b,c, ... , m, where all particles are of the same length, 
equal to the depth of the fuel bed, and are evenly spaced at distance 
L. Vf Na, No, Ne . . -, Nn are the number of particles of respective 


fuel types per unit area, and N is the total number of particles, then 
in a given length of bed each kind will be present in amount propor- 


' N. Ny N. N, 
tional to N? ry ‘aoe N’ If these fractions are represented as 


Rata tan + s » Se OT RR, ... . B repent the 


respective rates of spread through the different fuel types, the gross 
rate of spread of the mixture will be ' 


-. 
R=y—y cnet upa 


Rt mat +z (33) 


10 Equation (33) is analogous to one for total conductance of a series electrical 
; circuit of length L composed of several electrical conductors of length 1, each of 
different conductance but occurring n times. 


For example, let Na, Ns, Ne. . -» Nx be the number of times conductors of 
length 1, having conductances C,, C,, , C,, respectively, occur. The 
total length of the circuit will be L= N, a. Ni+N fe. + EN 
and the total conductance 

Ces * ot No+ nt Soe I 
Nz caries. 5 
c a+ Ott G f+. - + Cc. 
Then let 
'< SS eetome ne oa) be nas 
Pe Tthe .. +R RE 
where X,, Xn, X-, - - - » Xn are the fraction of the total conductance contributed 
by conductors type a, b,c, . . . , n respectively 
Then 
C= _— 
ee AR xs 


ee ae + 
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Let W., Wi, W., . . . » Wn be the total dry weight of the several 
particles per unit area of the fuel bed. Then 


W, , W, 

Agyad aaa? Oe a etc. (34) 
where A’s and y’s are the cross-sectional areas and the densities of the 
dry particles, respectively, and d is their length. 

| ae 


N=. 


By definition X= N? =H? ete. (35) 
and N=N,1+M+Ne+ ... +N, (36) 
Equation (35) becomes, after substituting (34) and (36): 
. Witwye 
<.=W. WwW. . W W, ete, (37) 


yet nag or 


2S OD rn. 


AnYn 


Substituting equation (37) in (33) and letting 
W, fy W, , W. We. 


= ‘ ae eae € 
‘ Aaa Arvo tact T Aste (8) 
the following equation results: 
i. 
R=, + oa ae (39) 
x, wv, i 
Agyaka Ary, ake Any, 


Equation (39) is an expression for the rate of spread of fire through a 
heterogeneous bed in terms of fuel composition and the rates of 
spread through each component fuel. The rate for each component 
is calculated by equation (32) with the value of the compactness, 
1/A, equal to that of the mixture. 

To calculate the rate of fire spread in a heterogeneous bed by 
equation (39) the spacing, L, which is not directly measurable, must 
be determined from the fuel and fuel-bed characteristics. The 
measurable factors from which Z may be computed are V,,V2, W, y, ¢, 
and \, where \ is determined by equation (15). 

For a given bed, the ratio of volume of voids to fuel-surface area is 
expressed by equation (15); the volume of a fuel bed by equation (16). 

The volume occupied by the fuel in the heterogeneous bed of volume 


V, is Va=V.t+V,t+V.+ --+ + Ve where V,, Vi,---, Ve are 
volumes of the different hinds: of eae It follows that 
W. 4 
vt ribo Sai oe (40) 


The total surface area es on ney in on bed is 


Wa, sWe , ocW, aaW, 
Ya + Yo T Yc . + Yn 


S= (41) 
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Substituting equations (16), (40), and (41), in (15) gives 


pay (Seige... 4 Me 
A=- a anes | Vee Ce ct te Sey Yn (42) 
oaWe Wa, ooW> , oWe oWe . «i San “sae 
Ya Yo c Yn 


Solving equation (42) for LZ’, then multiplying and dividing each term, 
the right side of the equation by the number of particles in that type 
of fuel, namely, Na, No, Ne +++ Nn, letting A,, A,, A., eo 


equal Wa , ny peers mc and assuming a length of 
vaNad yoN rd yNd YnNnd 





fuel particles equal to the depth of the fuel bed, the following equation 
results: 


L=[NMogAgXat opArXy+ocAcXe+ a ee + o,A,/ a) 
+ (AgXe+t AX, t+AXe+ eee +A,X,)]"" (43) 


where X,, X», Ac, + + +, Xn are defined by equation (37). 


By substituting the equations for X,, X,, X., + + -,.X, and simplify- 
ing, equation (43) becomes: 


iw [e+ $m oW, aos. east 4 onWn 


Yc 
, W,\ 7% 
Te m) (44) 








where B is defined by equation (38). 

Equation (44) is an expression for the particle spacing L for a 
heterogeneous fuel bed in terms of its composition. The assumption 
of fuel-particle length equal to depth of the fuel bed requires, of course, 
use of the coefficient C, in application of this measure to natural 
fuels. This coefficient is present in equation (32). 


EXPERIMENTS IN NATURAL HOMOGENEOUS FUEL BEDS 


EXPERIMENTAL PROCEDURE 


For evaluating the coefficient C,, ponderosa pine needles, being 
the most readily available uniform fuel in this region (California), 
were chosen as the natural fuel for the standard bed. In order to 
exercise control over all the variables, fires were studied on a model 
scale. Beds of needles 2 inches deep, 3 feet wide, and 8 or 12 feet 
long were prepared in trays with controlled compactness. These 
beds were placed in a wind tunnel (fig. 3) —— elsewhere (9) 
and burned on the horizontal, i. e., at zero slope. ure 4 shows the 
position of initial ignition of a fuel bed in the he Ignition was 
accomplished by an electrically heated coil placed on a match head 
inserted in the middle of a thin pyroxylin disk of 2-inch diameter. 
The fire traveling forward along the length of the bed was observed 
(fig. 5) from an overhead window. Each fire was allowed to burn 
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Figure 3.—Wind tunnel used in fire-behavior studies; test section is 30 feet long 
with a cross section 6 x 6 feet, over-all length 55 feet. 


Figure 4.—Position of initial ignition of a fuel bed in the tunnel. 
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until it reached approximately 18 inches in width and was then in- 
stantly suppressed by water. A typical burned area is shown in 
figure 6. 

The air velocity recorded for each fire was the velocity at 1 foot 
above the surface of the bed prior to ignition. A sample of fuel was 
collected at random-from the surface of the bed just prior to burning 
and its moisture content was determined by xylene distillation (4). 
The fuel temperature was taken to be equal to the ambient air tem- 
perature, measured when the moisture sample was taken. Approx- 





Figure 5.—Pine-needle fire 30 seconds Figure 6.—Typical burned area of pine 
after ignition, photographed from 
above; wind velocity 8 m. p. b. 


needles. 


imately 200 fires were burned in fuel beds of ponderosa pine 
needles (c=128 inches"') under varying conditions such as: wind 
velocity 2:0 to 12.0 m. p. h.; compactness, 1/\, 6.25 to 16.66 inches; 
moisture content 4.0 to 15.0 percent; air temperature 50° to 85° F. 
In order to simulate a standing type of fuel, such as brush, the beds 
were prepared by standing uniform twigs vertically and at equal spac- 
ing in sawdust treated with a fire-retardant chemical (fig. 7). The twigs 
were cut 7.5 inches in length from dead branches of ponderosa pine 
and segregated into three nominal size classes \, *{s, and 4 inch in 
diameter. The bed was set afire by a 30-gm. bundle of pine needles 
and to avoid the influence of the burning needles, the twigs were 
allowed to burn a distance of about 4 feet before zero time was called. 
679977—46—3 











110 


Journal of Agricultural Research Vol. 72, No. 3 


The fire was not extinguished until it had reached the end of the bed, 
which was 12 feet in length (fig. 8). Observations were made and 
data taken only in that portion of the burn beyond the first narrow 
region influenced by the mode of ignition. 


RESULTS 


Numerical values of C, were computed for each pine-needle fire by 
equation (32) using the observed rates of spread with C equal to unity. 





Figure 7.—Bed of ponderosa pine twigs ah er 
3%, inch in diameter with bundle of FIGURE a area of 
needles in the foreground used for Bs. 
ignition. 


The average values of C; were computed for several wind classes and 
are shown in figure 9. A value for the rate of spread of each fire was 
calculated by using appropriate values of C, from the curve in figure 9. 
Figure 10 shows the computed versus the actual rates of spread 
plotted to reveal the agreement. It will be noted that the points fall 
close to and evenly about a 45° line which represents agreement of 
calculated with experimental values. 

The experimental data for the twig fires are presented in table 1, 
together with the calculated values for the coefficient C. The values 
of C for individual fires vary considerably from its mean of 1.61 owing 
mainly to the experimental technique used in igniting the twigs and 
to the error in judging the exact position of the fire at zero time and 
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the time when the fire reached the end of the bed. Therefore, the 
value of 1.6 for C, for this type of bed, must be considered as only 
approximate until more data have been obtained. Since C was 
assumed to be unity for the pine-needle bed, the value of 1.6 is in the 
right order of magnitude for a fuel bed with twigs standing vertically, 
inasmuch as it was assumed that the flow of the gases passing the 
particles in the standard pine-needle bed was at right angles. It 
would have been more nearly correct, however, to assume that the 
gas flow is at right angles for only half the particles, and parallel or 
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WIND VELOCITY, V, (MILES PER HOUR) 


Figure 9.—Relation between temperature-ratio coefficient, C;, and wind velocity, 
V (miles per hour) at 1 foot above fuel surface, zero slope, and 1 minute from 
start of fire. Fuel: Ponderosa pine needles. 


longitudinal for the other half. Actually, in evaluating C, by the 
type of bed chosen as standard, it might have been more proper to 
assume a value less than unity for C, in which case the numerical 
values of C, would have been slightly higher. Then, with higher 
values of C, for beds with twigs vertically oriented where cross flow is 
more nearly approached, C would have more nearly approached unity. 
Practically, however, since C may vary with the technique used in 
sampling fuel moisture and measurement of wind velocity, it would 
probably differ from unity for other types of fuel beds. Consequently, 
for the sake of simplicity it would be better to assume cross flow for 
all types of beds and let C, evaluated in that way; be a constant 
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relating a particular kind of natural bed with a standard pine-needle 
bed, for which f. was calculated on the basis that the gases flow at 
right angles to the needles. 


TaBLE 1.—Experimental data and calculated values of C for model fires burning in 
fuel beds of ponderosa pine twigs placed vertically and spaced equally 






































Experimental data Calculations ! 
Fuel r . 
Fire | Wind | FUel| mois- | Fuel | Fuel | Rate 
(pum 10 |pera-| TWEE | olume| SBR | 
wei ware | tent | a be oL Seth yCp inf ¢ (=) | c 
| 
(V) | (to) | (M) (¢) (L) | (R) 
B. t. u./ B- t. u./ 
Inches Sq. ft. hr. Cu. ft. 
M p.h.| °F. | Pet. (2) |Inches| Ft./hr. °F. or. 
22 4 80 12.4 20.7 | 1.50 195 | 31.0 14.7 24.6 0.180 | 1.89 
3 4 86 9.1 20.7 | 1.50 223 | 31.0 14.7 22. 2 -176 | 1.91 
11 4 88 10.3 20.7 | 1.50 209 | 31.0 14.7 23. 1 | .174 | 1.84 
51 6 69 11.9 15.7 | 1.50 225 | 23.6 15.5 23.9 . 1389 | 2.04 
10 6 89 10.6 20.7 | 1.50 265 | 31.0 17.1 23. 4 . 125 | 1.46 
21 6 79 12.1 20.7 | 1.50 285 | 31.0 17.1 24.3 - 183 | 1.72 
1 6 87 9.8 20.7 | 1.50 335 | 31.0 17.1 22.7 . 126 | 1.81 
37 6 68 11.4 27.8 | 1.25 367 | 34.8 19.0 23.5 . 140 | 1.88 
49 8 68 12.5 15.7 | 1.75 222 | 27.5 17.0 24.4 - 105 | 1.22 
56 8 73 10.9 15.7 | 1.75 264 | 27.5 17.0 23. 3 -101 | 1.33 
65 8 62 11.3 15.7 | 1.63 257 | 25.5 17.0 23.3 - 109 | 1.51 
40 8 71 10.7 15.7 | 1.50 294 | 23.6 17.0 23. 1 .103 | 1.74 
52 8 68 12.2 15.7 | 1.50 242 | 23.6 17.0 24. 2 -105 | 1.53 
39 8 71 10.6 15.7 | 1.2 306 | 21.6 17.0 23.0 . 103 | 1.97 
48 8 61 11.4 15.7 | 1.25 244 | 19.6 17.0 23. 4 -110 | 1.88 
64 8 62 9.7 20.7 | 1.25 360 | 25.9 19.0 22. 1 -109 | 1.76 
47 8 59 11.9 20.7 | 1.13 285 | 23.3 19.0 23.7 -111 | 1.69 
2 8 87 9.5 20.7 | 1.50 492 | 31.0 19.0 22.5 091 | 1.71 
12 & 73 10.4 20.7 | 1.50 353 | 31.0 19.0 22.9 .101 | 1.39 
20 8 79 12.3 20.7 | 1.50 395 | 31.0 19.0 24.5 .097 | 1.59 
42 8 70 10.9 20.7 | 1.38 430 | 28.5 19.0 23. 2 . 103 | 1.90 
41 8 71 1L.1 20.7 | 1.25 339 | 25.9 19.0 23. 4 .103 | 1.66 
61 8 58 11.0 20.7 | 1.25 331 | 25.9 19.0 26. 4 .112 | 1.99 
43, 8 74 11.6 27.8 | 1.38 457 | 38.2 21.2 23.8 -101 | 1.36 
45 8 61 11.0 27.8 | 1.38 369 | 38.2 21.2 23.1 -110 | 1.16 
63 8 62 8.7 27.8 | 1.38 497 | 38.2 21.2 21.3 . 109 | 1.42 
44 8 74 10.2 27.8 | 1.25 613 | 34.8 21.2 22.8 | -101 | 1.91 
46 8 60 11.4 27.8 | 1.13 497 | 31.3 21.2 23.3 | -110 | 1.92 
50 8 67 13.3 27.8 | 1.00 395 | 27.8 21.2 25.0 | - 106 | 1.78 
57 8 73 11.0 27.8 | 1.00 465 | 27.8 21.2 23.3 | -101 | 1.86 
lla 5.8 71 11.6 20.7 | 1.38 222 | 28.5 16.7 23.8 | . 141 | 1.57 
12a 5.8 71 11.7 20.7 | 1.25 233 | 25.9 16.7 23.8 | - 141 | 1.81 
13a 5.8 71 11.0 15.7 | 1.25 171 | 19.6 15.1 23 3 | . 141 | 1.90 
l4a 5.8 71 11.1 15.7 | 1.50 149 | 23.6 15.1 23. 4 | . 141 | 1.38 
18a 5.8 8&8 4.6 27.8 | 1.25 457 | 34.8 18.7 18.6 | . 130 | 1.70 
19a 5.8 89 5.3 27.8 | 1.13 436 | 31.3 18.7 19. 2 - 129 | 1.84 
20a 5.8 89 5.5 20.7 | 1.50 320 | 31.0 16.7 19. 4 . 129 | 1.55 
25a 5.8 94 | 5.0 15.7 | 1.75 236 | 27.5 15.1 19. 1 - 125 | 1.36 
26a 5.8 92 4.6 20.7 | 1.13 297 | 23.3 16.7 18.7 .127 | 1.81 
33a 5.8 91 5.7 27.8 | 1.38 395 | 38.2 18.7 19.5 } -127 | 1.37 
34a 5.8 89 5.4 15.7 | 1.38 222 | 21.6 15.1 19.3 | - 129 | 1.69 
35a 5.8 89 5.0 15.7 | 1.63 242 | 25.5 15.1 19.0 - 129 | 1.54 
40a 5.8 84 | 7.8 | 1.00 347 | 27.8 18.7 20.6 | . 132 | 1.82 
59a 5.8 81 7.2 20.7 | 1.38 77 | 28.5 16.7 20.6 | . 134 | 1.61 
27a 7.8 92 5.1 15.7 | 1.75 391 | 27.5 16.7 19. 1 | O91 | 1.46 
28a 7.8 91 5.4 20.7 | 1.38 395 | 28.5 18.7 19.3 | . 092 | 1.32 
4la 7.8 84 6.4 15.7 | 1.38 250 | 21.6 16.7 20.0 | .097 | 1.34 
45a 7.8 85 6.6 15.7 | 1.63 244 | 25.5 16.7 20. 2 | .096 | 1.12 
46a | 7.8 85 7.0 15.7 | 1.50 256 | 23.6 16.7 20. 5 | .096 | 1.28 
| | 





| 
| 
| 
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nsCanf c (2 =) ] 
1 According to equation (32) C=——— f""i7 4 where fuel density, y, =26.1 Ibs./eu. ft.; tr= 


1,500° F.; t:=550° F. 
2 1+number of inches, 
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Figure 10.—Calculated versus experimental rates of forward spread for ponderosa 
pine-needle fires. 


VARIABLES MEASURABLE IN FIELD PRACTICE 


Equation (32) has been derived for the rate of spread in terms of 
the fundamental variables, which in practice are generally wind 
velocity, moisture content, time, and slope. Fuel-type or cover-type 
classifications also are used almost universally in recognizing the 
rate-of-spread differences attributable to fuel size and cover density, 
which the present study has revealed to be important and significant. 
Compactness, which incorporates fineness and crowdedness of the 
particles in the fuel bed, is used in this analysis to define cover density. 
Two other important measurable variables have been added, namely, 
fuel temperature and fuel-particle density, heretofore not generally 
recognized or applied in practice. 

These measurable variables influence fire spread only as they pro- 
duce changes in the fundamental variables. Some operate on as 
many as three or four fundamental variables. For example, the 
fuel temperature, f,, affects film conductance for convection and con- 
duction, f,., heat-transfer factor, f,, specific heat of moist fuel, ©,, 
and in addition is itself a fundamental variable in the spread equation. 

Figure 11 shows curves for rate of spread as a function of all the 
variables, except slope and time, as interpreted through application 
of the fundamental relationships expressed in equation (32). 
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Figure 11.—Effect of measurable variables on rate of forward spread of fires. 
For each graph five variables are held constant while the sixth is varied over a 
practical range. When held constant, the variables have values as follows: 
V-6.0 m. p. h.; 1/A=10 inches -!; y=31.6 lbs./cu. ft.; 4=100° F.; c=128 
inches ~!; M=8 percent. 
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INFLUENCE OF MEASURABLE VARIABLES 
EFFrect oF WIND VELOCITY 


Wind velocity influences fire spread by causing changes in the 
magnitude of film conductance, f., and temperature ratio coefficient, 
C,. The film conductance varies approximately with the square root 
of the wind velocity (see equation 12). Since rate of spread is directly 
proportional to f., it will vary with the square root of the wind ve- 
locity. Most of the increase in fire spread that is due to the wind 
velocity, however, results from its effect on the temperature ratio 
coefficient, C,. Increasing velocity decreases the flame angle, and 
this results in more effective preheating of the fuel; that is, fuel 
temperature, ¢t,, increases with wind velocity. From equation (29) 
it will be seen that higher temperature values ¢, of the fuel particle 
about to ignite will give smaller values of C,; thus C, is some inverse 
function of wind velocity. The way C, varies with wind velocity is 
shown in figure 9. Since the rate of spread varies inversely as the 
logarithm of C,, it will vary directly as the logarithm of the velocity. 
With wind measured at 1 foot above the fuel bed, the product of 
these separate functions indicates that rate of spread is approximately 
proportional to the first power of wind velocity for velocities less 
than 5 m. p. h. and to the 1.5 power of wind velocity for velocities 
from 5 to at least 12m. p.h. The upper velocity limit (to which the 
1.5 power holds) has not been definitely fixed. The first power for 
wind velocities less than 5.0 m. p. h. is confirmed by experimental 
results obtained by Curry and Fons (6) from test fires under field 
conditions. 


EFFECT OF FUEL S1zE 


In this analysis, the surface-volume ratio, o, of the fuel particle 
has been used to express the fuel size. The presence of o in the numer- 
ator in equation (32) does not necessarily imply that rate of spread is 
directly proportional to c. In fuel beds of constant compactness, 
variations in o will affect the distance Z between fuel particles (equa- 
tion 21). In addition, o affects the film conductance (equation 12). 
The effect of o on ignition temperature, ¢;, has not been considered in 
the analysis, but significant variations of ¢; with o are expected 
(equation 10). By holding all other factors constant, the expression 
for rate of spread in terms of o becomes: 


R= (ao+bo°+-¢) (1+do)’* 


where a, b, c, and d are the proper constants for the given values of 
wind velocity, initial fuel temperature, compactness, and moisture 
content. For fine fuels this expression very closely approximates a 
straight line having positive slope. Experimental results (7) show 
that the rate of spread for fine fuels is very nearly inversely propor- 
tional to the fuel size, i. e., it varies directly as the surface-volume 
ratio of the fuel particle. 


EFFECT OF COMPACTNESS 


For want of a more descriptive word, the variable incorporating 
fineness and crowdedness of the fuel particles is here termed compact- 
ness. The reciprocal of the ration of volume of voids to the surface 
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area of the fuel in the bed, denoted by \, may be used as a direct 
measure of this element. By definition, the spacing of fuel particles, 
L, varies as the square root of X. Since rate of spread is proportional 
to L, it too will vary as the square root of \. For a fuel particle of 
given size it is necessary to increase the distance L between particles 
to obtain greater values of (equation 21); this would indicate that 
a unit increase at large values of \} would not increase the rate of 
spread as much as at the small values. There must be, of course, a 
limiting value of \ at which the fires would fail. Beyond a critical 
limiting point further separation of fuel particles decreases the volume 
of flame, and the hot gases produced in diminishing quantity are 
cooled and dissipated to such an extent that they are no longer 
effective in igniting the unburned fuel. Experimental results (7) 
have confirmed the square root of \ relationship. 


EFFEcT OF FUEL TEMPERATURE 


Changes in the initial temperature of the fuel, tf), affect the four 
fundamental variables, film conductance, heat-transfer coefficient, 
specific heat of the moist fuel, and temperature ratio v: e The 

tS 3 

film conductance changes with the fuel temperature because the kine- 
matic viscosity and thermal conductivity of the film are affected 
(equation 12). An increase of 50° F. in fuel temperature, tj, will, how- 
ever, decrease the film conductance only 0.5 percent and may be dis- 
regarded. The heat transfer coefficient, f,, is also negligible for large 
changes in fuel temperature, because flame temperature is large in 
comparison with ¢t, (equation 14). The specific heat of the moist fuel, 
C,, increases with a rise in fuel temperature (equation 27). Since 
rate of spread is inversely proportional to C, (equation 32), an increase 
in C, decreases the rate of spread. The important effect of fuel tem- 
perature on fire spread is that a change in temperature changes the 
ignition time, 6;, as shown by equation (10), with 


—ft . wat! . 
aes ) replacing wae - Changes in fuel temperature affect 
liek ' me S| 


the magnitude of the ratio . 
mae | 


produces a comparatively large change in In | o(? a =) The rela- 
tive extent of the fuel temperature’s effect on fire (rr may be deter- 


mined by calculating the magnitude of the factors 1/ln [ce "| 
mer y 


» and a small change in this ratio 


Set, - P 5 ae sie 

and = in equation (32). For conditions indicated in figure 11, an 
YUp 

increase in et temperature from 70° to 120° F. causes the factor 


t 
1/ln [eGe _ » and consequently rate of spread, to increase 35 per- 


; +f, ei 
cent; with en same change, the factor 1222 " causes a decrease of 5.5 


ee 
percent. The total effect of the 50° F. increase in initial fuel tempera- 
ture, under the conditions described, is thus a net increase of 29.5 per- 
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cent in fire spread. Since C; decreases with wind velocity, a larger 
increase in spread results with rising wind velocity for a given tem- 
perature change; i. e., for the same conditions as above except that 
with wind velocity at 10 m. p. h. a rise in temperature from 70° to 120° 
F. produces a 94-percent increase in the rate of spread. The tendency 
of a fire to slow down soon after a shadow from a cloud is cast on the 
fuel ahead of the fire has frequently been observed. The fuel tempera- 
ture in the rate-of-spread equation offers a possible explanation of this 
phenomenon, since fine fuels in complete shade soon assume air-tem- 
perature conditions, while in direct sunlight a much higher fuel 
temperature generally prevails. 


EFFECT OF FuEL DENsiTy 


The density of the moist fuel, y, appears in the form of an inverse 
proportion as one of the fundamental variables in the expression for 
rate of spread. Since the density of a body is its mass per unit volume, 
more heat is required to bring a body of higher density than one of 
lower density to a given temperature. Consequently, the ignition 
time for fuel particles of a given volume is directly proportional to 
density (equation 10), and by the same token fire spread is inversely 
proportional to density. 


Errect or .MotisturE CONTENT 
M 

1+M 
in the equations for density and specific heat of the moist fuel, respec- 
tively. The factor 1+ in the two expressions cancels, leaving M in 
the product yC,, which is an expression for specific heat based on 
volume or volumetric specific heat. Since fire spread varies inversely 
with yC, (equation 32), it will also vary with the fuel-moisture con- 
tent. Moisture content can be expected to have the same effect on 
the rate of spread as density has, in that the greater the moisture the 
more heat is required to bring the fuel to the ignition point and the 
slower it will burn. Since rate of spread must become zero at some 
given moisture content for a particular fuel size, it is doubtful whether 
the effect of absorbed water on rate of spread is accomplished entirely 
through its effect on y and C, for the complete range of moisture 
content. At high moisture contents it would seem probable that the 
excessive water vapor acts as a diluent for the inflammable fuel gases 
to the extent that rate of spread decreases more than y and C, might 
indicate. If the moisture is sufficiently high, the water vapor may 
dilute the gases to the point where combustion is entirely prevented; 
this would determine the upper moisture content limit. For the range 
of moisture content studied (4 to 15 percent) equation (32) is valid; 
otherwise values of C, computed for the experimental fires for any 
given wind velocities would have varied with the moisture. This, 
however, was not found; thus the effect of absorbed water on rate of 
spread is apparently accounted for by y and C, for moisture contents 
up to 15 percent. 


The moisture content, 17, appears in the factors 1+.M and 


EFFECT OF SLOPE 


Because hot gases or convection columns tend to rise vertically 
with respect to a horizontal surface, an increase in slope of terrain 
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decreases the angle between the fuel bed and hot gases. Slope, there- 
fore, has the same effect as wind in bringing the flame closer to the 
fuel. As such, it may be considered an added component of wind 
velocity. As the flame angle decreases, the fuel starts to absorb 
heat, both by radiation and convection, farther ahead of the actual 
burning front than with a flame perpendicular to the fuel surface. 
This results in increasing ¢,, temperature of the fuel particle about to 
ignite, which in turn makes C, smaller (equation 29). Since the rate 
of spread varies inversely as the logarithm of C,, it will vary directly 
as the logarithm of the 

TIME = 0 slope. The effect of 
iia slope on rate of 

\ spread indicated by this 
ba relationship is shown in 
' figure 12, where the 

~ lines drawn for differ- 
a n,n ee ent wind veloci- 
ie >—~,4—_ a — = — ties are broken to indi- 


—_- 
_— 

— —. 
—_— —. 


~~ hee = === === cate that their true 

3 —— — positions have not yet 
been experimentally 
established. 


TEMPERATURE-—RATIO COEFFICIENT, © 


EFFECT OF TIME 


SLOPE, J, (PERCENT) Immediately after 
the start of a fire the 
volume of hot gases at 
its head increases pro- 

pias a" gressively. As the vol- 
ume of the flame rises 
the rate of heat trans- 








we fer by radiation to the 
AON fuel particles is in- 
\ i aT i I ° 
\ iy, Ze creased ; t h us t, in- 
eg a Se eran seks creases with time until 
SN > +477: the volume of hot gases 
ne eeran ener emeerecuces aac ce from the flame front 
een ettancigcnvectea has reached a more or 
less constant value. 
As in the case of slope, 
TIME, 8, (HOURS) increasing values of t 
Figure 12.—Curves showing how the tempera- result in decreasing 
ture-ratio coefficient may be expected to vary Values of C, (equation 
with slope and time for various wind velocities, V. 29); the rate of spread 
will therefore vary as 
the logarithm of time. This relationship was also indicated by earlier 
studies of field test fires (6). Model fires are of too short duration to 
permit studying the effect of time on C;. Data from which to derive 
actual values of C; at different time intervals must therefore be ob- 
tained from larger scale field experiments, which have not yet been 
undertaken. The way in which C, is expected to vary with time for 
several wind velocities and given slope is illustrated in figure 12. 
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VALUES FOR GROUPS OF FUNDAMENTAL VARIABLES IN TERMS OF MEASURABLE 
VARIABLES 


In table 2 are presented values of ¢, y, and EF for different kinds of 
fuel. In tables 3 to 6, inclusive, are computed values of groups of 
fundamental variables that appear in equation (32), as functions of 
those measurable variables affecting the group. 


TABLE 2.—Surface-volume ratio, perimeter, density, and — factor of several fuels 























Surface- cg : : 
volume — Density | <a 
Fuel type aie 
| o yg | Y | E 
Ponderosa pine: 1Inches-!| Inches | Pounds/cu. ft. 
TH AN... ewan = ce none se en ore : J 4/D 2eD 29. 1 3.6 
eee aS pe Pe 4/t 3 4t 26.7 4.0 
Needles: 
EEE EEE LALOR eee ELIE 128 131 31.6 4.1 
ag ieenigdinmanewsc ance 165 . 108 33.5 | 4.2 
NS Ree eee nae, La ese eE I 140 . 136 | 35. 4 | 4.3 
_ (ee hy eee eee eee 7 nskekel 124 . 165 | 34, 3 4.5 











1 1+number of inches. 
2 Diameter of fuel. 
3 Thickness of fuel. 


TABLE 3.—Values! of f-+f,, (B. t. u. si. ft. hr. °F.) as a function of wind velocity 
and surface-volume ratio of fuel 





: | Wind velocity, in miles per hour indicated at 1 foot above fuel bed 
Surface-volume ratio | 
(inches ~!)2 





| 





| 2 4 6 8 10 12 
RON | 9.7 | 11.8 | 13.5 | 14.8 | 16.1 | 17.2 
20 | 11.5 | 14.5 | 16.8 | 18.9 | 20.7 | 22.2 
ESS Sie Bee | 14.3 18.5 | 21.8 | 24.5 | 26.9 | 29.1 
ee Seton es | 17.6 22.6 | 26. 8 30.7 | 33, 8 | 36. 8 
“TAOS 2 is Se ni aa 20. 3 26. 3 | 31.3 35.2 | 39.0 | 42.3 
(aaa 22.7 29.4 | 34.9 39. 2 | 43.6 | 47.2 
160 25.1 32. 2 | 38. 4 | 43.5 | 47.9 | 51.8 








1 Values in table are within +0.5 percent for a range of fuel temperature from 50° to 140° F. 
2 1+number of inches. 


" oL * a 
TABLE 4.— Values ! of Ewa function of surface-volume ratio of fuel and fuel bed 
4 


compactness 





Fuel bed compactness (inches -!)? 

















Surface-volume ratio tad 
(inches ~!)2 | | | 
iit Bch CARD tee, ee We ee a a ee 
| | | 
7.14 4, 58 3, 32 245| 1.87| 1.63| 1.41| 1.22 
10. 05 6. 40 4. 58 3. 32 | 2. 45 | 2.08 | 1.73 | 1.41 
14, 18 9. 00 6. 40 | 4. 58 3. 32 2.77 | 2. 24 | 1.73 
18, 7: 11. 87 8. 43 6.00 | 4. 30 3. 56 | 2. 83 | 2.12 
22. 38 14, 18 10. 05 7.14 5.10 | 4.21 3. 32 | 2.45 
25. 51 16. 16 1.45] 812] 5.79] 4.77 3.74 | 2.74 
28. 30 17.92 | 12. 69 | 9.00 | 6. 40 | 5. 26 | 4.12 3.00 





- ; : L ; 
1 To obtain oL for equation (32), multiply by an appropriate shape factor, E, of the fuel. 
2 1+number of inches. 
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a 3 P . : 
TABLE 5.—Values! of 1? as a function of moisture content and temperature of fuel 
7 Pp 4 





Fuel moisture content (percent) 





Fuel temperature | 
(°F.) 











| 2 | 4 | 6 | 8 | 0 | 12 “ | 16 
| | | 
eek ee 0. 546 0. 597 0. 648 0. 699 0.750 | 0. 801 0. 852 0.903 
Ui re eas . 556 . 608 . 660 .713 - 765 | .817 .870 “922 
“Aaa Oe ewe rmeaaied gmat 620 | . 674 727 . 781 | 835 | 888 * 942 
ee eee | .577 682 2 | . 687 . 743 | .798 | 853 . 908 963 
Si eaaaly _ 588 : . 702 . 759 . 816 872 929 986 
eae aia | 600 | ‘6 4 717 .776 . 835 . 893 952 1.010 
1 








os * 4 Cp ° 
1 To obtain yCp for equation (32), multiply — by the density of the fuel. 


P - tr—t, e ; A A 
TaBLE 6.—Values'! of ln [c. C= ) | or surface fires as a function of wind velocity 


and fuel temperature, to, for time 1 minute from start of fire and zero slope 





Wind velocity (in miles per hour indicated) at 1 foot above fuel bed 
Fuel temperature (°F.) 











| 2 | 4 | 6 | 8 | 10 12 
40 ‘iin edksVetaivebedGased 0. 2904 0. 2078 0. 1596 0. 1249 0. 0971 0. 0760 
| BERS SRS eee Rees a cere cna . 2769 . 1939 - 1458 - 1106 . 0834 - 0620 
Dt ininatblictebeass<tanwen - 2631 . 1798 . 1319 . 0971 . 0695 . 0478 
100 Scape ehmeaees . 2484 . 1664 . 1178 . 0825 . 0554 . 0344 
120 ibcumccndasnanebneee . 2343 . 1519 . 1035 . 0686 . 0411 - 0198 
a Oe thawakdvahwéasit . 2199 1371 . 0889 0535 - 0266 . 0070 

| 





1 Flame temperature, fy, and ignition temperature, f;, are considered constant at 1500° and 550° F., re- 
spectively. 


SUMMARY AND CONCLUSIONS 


The analysis of rate of fire spread in light forest fuels reported in 
this paper is based on the theory that the spread of a fire can be ex- 
pressed as successive ignitions of adjacent fuel particles and that its 
rate is therefore governed by the time required to raise successive 
fuel particles to ignition temperature. Equations have been derived 
for rate of fire spread in homogeneous and heterogeneous fuel beds, 
taking‘into account the physical characteristics of the fuel particles, 
the arrangement of the bed, and the pertinent attributes of the atmos- 
phere. The spread equations are expressed in terms of the following 
fundamental variables: (1) Film conductance; (2) heat transfer factor 
for radiation; (3) ignition temperature; (4) fuel particle spacing; (5) 
surface volume ratio of fuel; (6) specific heat; (7) density of fuel; and 
(8) fuel temperature. It has been shown that the variables ordinarily 
recognized in field practice—wind velocity, fuel moisture, fuel density, 
fuel size, fuel bed compactness, slope, and time—influence fire spread 
only as they produce changes in the fundamental variables. 

The results of field and laboratory experiments were used to check 
the theory and analysis, as well as to provide quantitative information 
for determining the magnitude of the effects of variables measurable 
in field practice on the rate of spread. 
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A COMPARISON OF THE CHEMICAL COMPOSITION OF 
DIPLOID AND TETRAPLOID CORN! 


By G. H. E.uis, research chemisi, United States Plant, Soil, and Nutrition Labora- 
tory, LL. F. RaNpouwrn, cytologist, Division of Cereal Crops and Diseases, Bureau 
of Plant Industry, Soils. and Agricultural Engineering, and G. Marrone, 
chemist, United States Plan, Soil, and Nutrition Laboratory, Agricultural 
Research Administration, United States Department of Agriculture 


INTRODUCTION 


In crop production studies, increasing attention is being given to the 
evaluation of the crops produced in terms of constituents of nutri- 
tional significance, as well as in terms of total yield. Although exten- 
sive animal feeding trials are necessary to evaluate a crop completely, 
chemical data are useful indicators of the nutritional value of plants. 

In 1940 Randolph and Hand (10) * reported that doubling the num- 
ber of chromosomes in pure yellow corn resulted in an increase of 
approximately 40 percent in the total carotenoid pigment content and 
in the biologically active provitamin A fraction of the grain. Further 
analyses on both the grain and stover of similar strains of corn indicate 
that doubling the number of chromosomes may bring about other 
changes in composition that appear also to be of nutritional signifi- 
cance. 

MATERIALS AND METHODS 


Analyses of the grain involved comparable diploid and tetraploid 
cultures in each of four stocks of corn (Zea mays L.), as follows: 


A. A commercial inbred line of yellow dent corn, designated 4—8d, and an auto- 
tetraploid line derived from it through direct chomosome doubling induced by 
the heat-treatment technique of Randolph (9). 

B. A commercial yellow dent line of tetraploid corn that was inbred for seven 
generations after it was produced by the heat treatment of a diploid hybrid com- 
bining inbred lines of Webber Dent, Illinois A, and Luces Favorite, and diploid 
cultures originating from parthenogenetic diploid plants that arose spontaneously 
in the tetraploid line. 

C. A first generation hybrid of the tetraploid lines of stocks A and B. 

D. A white flint line of tetraploid corn that was inbred for four generations 
after it was produced by heat treatment of a diploid hybrid combining two similar 
lines of white flint corn, and diploid cultures originating from parthenogenetic 
diploid plants that arose spontaneously in this tetraploid line. 


Analysis of the earfree stover involved A, B, and D utilized in the 
analyses of the grain and the following four additional stocks: 


EK. A commercial yellow dent tetraploid line, inbred two generations after being 
produced from a synthetic diploid hybrid involving inbred lines of Webber Dent, 


1 Received for publication August 16, 1944. Cooperative investigation of the U.S. Plant, Soil, and 
Nutrition Laboratory, the Division of Cereal Crops and Diseases, Bureau of Plant Industry, Soils, and 
Agricultural Engineering, Agricultural Research Administration, U. 8. Department of Agriculture, and 
the Department of Botany, New York State College of Agriculture, Cornell University, Ithaca, N. Y 

2 Also professor of botany, New York State College of Agriculture, Cornell University. 

3 Italic numbers in parentheses refer to Literature Cited, p. 130. 
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Illinois A, Luces Favorite, and white flint, and a diploid culture that orignated 
from a parthenogenetic diploid plant occurring in this tetraploid line. 


F. Another tetraploid line of the same origin as E, inbred three generations, 
and a derived parthenogenetic diploid culture. 


G. A third line of the same origin as E and F, inbred three generations, and a 
derived parthenogenetic diploid culture. 


H. A tetraploid line of the same origin as B, inbred three generations and 
parthenogenetic diploid cultures derived from it. 

The above stocks A to H were analyzed in 1940 (tables 1 to 3). 
In 1941 analyses of the stover in stocks A, B, and D were repeated in 
order to determine the extent of seasonal variation with respect to 
chemical composition, the same stocks being given different culture 
numbers in each of the 2 years as shown in table 3. 

The samples of grain and stover were harvested at maturity, the 
grain being air-dried immediately after harvest to a moisture content 
of approximately 14 percent. The stover samples were taken from 
fully matured plants before they had been injured by frost or desicca- 
tion and while still green. After the stover samples had been air- 
dried, they were ground in a hammer mill. Suitable portions of this 
material and of the corn grain were prepared for analyses by grinding 
in a Wiley mill to pass through a 1-mm. screen. 

All of the results given are expressed as percent on a dry matter 
basis. The dry matter was obtained by drying in a vacuum oven at 
95° to 100° C. for 5 hours. 

Total nitrogen was determined by the boric acid modification of the 
Kjeldahl method. 

Ether extract, ash, and crude fiber were determined by the official 
methods of the Association of Official Agricultural Chemists (1). 

The cellulose determinations were made by the Crampton-Maynard 
method (3) and the lignin determinations by a 72-percent sulfuric 
acid method, the detail of which have not yet been published. It 
consists essentially in extracting the air-dried sample with an alcohol- 
benzene mixture, followed by digestion with pepsin to render protein 
material soluble. After additional acid hydrolysis, the residue is 
treated with 72-percent sulfuric acid. Water is added to bring the 
sulfuric acid concentration to 3 percent and boiled. The residue after 
filtration and correction for its ash content is considered to be lignin. 


RESULTS 


The values for the nitrogen content of the corn grain are presented 
in table 1. This table and succeeding tables show that more than one 
culture of each stock ordinarily was analyzed. These cultures within 
each of the stocks originated from seed of selfed ears which were 
produced by sister plants in the preceding generation. For example, 
in stock B cultures 612 and 613 originated from the selfed seed of 
spontaneous parthenogenetic diploid plants, and cultures 614 to 618 
from tetraploid sister plants of the same culture in the preceding 
generation. 

In all comparisons, the nitrogen content was higher in the tetraploid 
stocks than in the comparable diploid stocks. The largest increase 
was in stock D, a white flint corn. The mean value with its standard 
error for all the diploid values was 1.799 + 0.0887, while the mean value 
for the tetraploids was 2.087 + 0.0349. The difference of 0.288 + 0.0953 
is statistically significant at odds of 100:1 and represents an average 
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increase of 15 percent. In this comparison, the values of stock C are 
omitted. These values are of interest since they show that this 
hybrid of the tetraploid stocks A and B was considerably lower in 
nitrogen content than either of the parent stocks. 

The individual values recorded for stock C represent separate 
analyses and were included to indicate the adequacy of the sampling. 
The determinations were made on five replicate 20-gm. samples and 
show relatively close agreement among themselves. All the other 
samples were 50 gm. in weight with the exception of the tetraploids 
in stock A, which were 20 gm. 


TABLE 1 ore content slhallva corn seeing 1940 





ve | Culture | Diploid | erage | Culture || Tetraploid | 4. ,,, 
Stock No. Diploid | Average No.! | Tetra | Average 
| | 
| | | | 
| Percent | Percent | | Percent Percent 
Nee Eee oe 091 | Yas 562 || 2. 309 | 
| 565A 2,173 | 
2. 139 | 565B 2. 216 2. 233 
JR ea ee ee 612 2.016 614 2. 042 
613 1. 981 615 || 2. 224 
616 2. 042 
1. 999 618 || 2. 216 2.131 
RAE Ee oe, Sa GE CPR eens: ERA ate 875-1 || 1.714 
575-2 || 1, 758 
575-3 | 1.781 
575-4 1.778 
575-5 | 1, 728 1. 752 
OS cee oa et eee €01 | 1, 522 | 603A 2.042 
€02 | el a 603B || 2. 005 
€06 | Ce cae 604A | 2. 082 
607 LO foc. 604B || 2. 107 
608 | i" Deir 605 1. 960 
1. 651 609 2. 026 
610 || 1,931 


611 1.979 | 2.017 





! A and B in this column represent abies of the same culture. 


There was no significant difference between the tetraploid and diploid 
stocks with respect to ether extract, ash, or crude fiber, as indicated 
by the data presented in table 2. 


TABLE 2.—Composition of corn grain, 1940 





| 
| | 
Stock Culture 











— | 
. Ether Crude 
No.! Type extract Ash | Nitrogen fiber 
ory g 
Percent | Percent Percent Percent 
ee Meg. no. as cnn ene aon 3.96 | 2. 00 2. 139 2. 66 
i icacacecs tl 565A | | ace eee 3.75 | 1.90 2.178 3. 20 
| 565B ..do Blix dine sioheeiec ae keith 3. 60 2.00 2. 216 2.96 
| 601 | Diploid _- 5. 37 A ef 1. 522 1. 68 
} 602 | RS eee re cera 4.77 1.61 1. 523 1.77 
D | 603A | 7. etraploid_- RSS ae 5. 33 1. 66 2. 042 1.74 
Perea 603B |. 5 LIN SLID MPs, 5.40 1.59 2. 003 1.79 
|} 604A Lo GPa beta 5. 58 1. 60 2. 082 1.89 
| "2p Et eae eee 5. 56 | 1.70 | 2. 107 2.15 











' A and B in this column represent replicates of the same culture. 


The values obtained for the nitrogen content of corn stover for both 
1940 and 1941 are shown in table 3. The difference between the 
diploids and tetraploids with respect to nitrogen content was greater 
for the stover than it was for the grain. In all comparisons, with the 
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exception of stock G, the tetraploid stover was higher in nitrogen than 
the diploid stover. Thus, the mean value was 1.224+0.0372 for the 
diploid and 1.464 +0.0319 for the tetraploid in 1940, and 1.452 +0.0882 
for the diploid and 1.941+0.0272 for the tetraploid in 1941. These 
values represent an increase in the nitrogen content of the tetraploid 
stover as compared with the diploid stover of 20 percent in 1940 and 
34 percent in 1941. The seasonal difference with respect to the diploid 
and tetraploid cultures was greater for certain cultures than for others, 
but there was no apparent explanation for these differences. 


TaBLE 3.— Nitrogen content of corn stover 


















































rT 
1940 | 1941 
Stock | ® ' | | Stock On | % ' 
Culture|,);.);4|Culture| Tetra- | ‘ulture|7;);.):4|Culture} Tetra- 
No. ey No. | ploid || _ No. an No. ploid 
| | 
Percent | | Percent | Percent Percent 
' (ae 091 1,317 562 1.640 || A-- suse) 1. 445 66A 1. 699 
565 1. 573 || | 68B 1.746 66B 1. 826 
—— || 663A 1.742 
Average_- ‘ 1. 607 || 663B 1. 861 
| 664A 1.770 
ae 612 1. 296 614 1. 379 |} 664B 1.818 
613 1. 042 615 1, 365 ] 669A 1.747 
616 1, 299 | 669B 1. 880 
618 | 1.397 || 687A | 1.971 
Lesine unseen At 87B | 1.878 
Awetege..>4.--:.-.. 1.169 — 1. 360 || 688A 2. 000 
\| 688B | 2.077 
Se 601 1.117 603 1. 627 || | — Seana 
602 | 1.211 604 | 1.584 || = Average____. | 1.594 1. 856 
606 1. 283 605 1. 646 {| | 
607 1.173 609 1.238 || B . 409 | 1.458 407 1. 792 
608 1. 045 610 1. 230 1] 413 1.773 408 1. 661 
611 1.328 || 459 1. 366 410 1. 885 
— —— || 464 | 1.595 411 1.970 
Average___- 1. 166 oda 1. 442 |} 412 1. 962 
| 414 2. 224 
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Since the ears of the tetraploid plants were not as well filled as those 
of the diploid plants, a test was made to determine whether this fact 
had influenced the results. In the eight stover cultures designated by 
A and B (table 3, columns 7 and 9), samples were harvested from 
plants that produced no seed, the ear shoots having been protected 
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with glassine bags during the pollinating season, and from plants that 
did produce seed. The fertile plants are designated A; the sterile 
ones B. The ears in all cases were removed before analysis. The 
mean values with the standard error were 1.773 +0.0608 for the fertile 
and 1.914+0.0559 for the sterile corn. This difference in nitrogen 
content is not statistically significant. 

Certain samples of the corn stover that showed marked differences 
in the nitrogen content of the diploid and tetraploid cultures were 
subjected to additional analyses for ether extract, ash, crude fiber, 
lignin, and cellulose. The data are shown in table 4. The values for 
nitrogen content given in table 3 for these samples are included in 
table 4 for comparison but are expressed as crude protein. The factor 
6.25 was used to convert nitrogen to crude protein. The results are 
expressed as crude protein so that the relationship between cellulose 
content and the nitrogenous constituents of the stover can be more 
readily seen. The only values other than those for crude protein con- 
tent that show any appreciable differences are those for crude fiber 
and cellulose. Thus, in stocks A and D, the increased crude protein 
content is accompanied by a roughly corresponding decreased cellulose 
content. Where the increased crude protein is not balanced closely 
by a decrease in cellulose, a change in the amounts of other carbo- 
hydrates such as starch or the pentosans must occur. That the crude 
fiber determination is essentially a measure of the cellulose content 
of corn stover is indicated by the close parallelism between the crude 
fiber and cellulose values. Seasonal differences in cellulose and lignin 
content were not as pronounced as those for crude protein content. 


TABLE 4.— Composition of corn stover 
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DISCUSSION 


The data presented show clearly that tetraploid corn is higher in 
nitrogen and lower in cellulose than comparable diploid corn, and that 
both types have essentially the same amount of the fatlike constitu- 
ents included in the ether extract and the inorganic constituents of 
the ash. But it is possible that these results do not accurately reflect 
the actual effect of chromosome doubling on the chemical composi- 
tion of corn because of the influence of factors other than the differ- 
ence in chromosome number. However, the experiments were 
carefully designed to minimize genetic and environmental differences 
that might affect the results. Various stocks were selected for study 
in order to sample different types of corn and also to include a sufh- 
cient number of comparisons to permit average differences between 
the diploid and the tetraploid to be obtained that would be more 
significant than differences based solely on one or a few comparisons. 

Since stage of maturity can influence composition to a marked 
degree, the samples were harvested only from fully mature, green 
plants that had been grown in randomized plots in the same field. 
It is also possible that the tetraploid plants were more leafy than the 
diploids. If they were, this might have accentuated the observed 
differences, since the leaves probably contain more nitrogen per unit 
of dry weight than do the stems. Unfortunately, this ratio was not 
determined, but the diploid and tetraploid plants were not noticeably 
different in this respect. 

The transformation from the diploid to the tetraploid state resulted 
in an appreciable change in the direction of higher nutritive values 
so far as can be judged by chemical analysis. Since the question of 
the significance of chemical analysis in evaluating feedstuffs has been 
adequately discussed elsewhere (3, 8), it is sufficient to point out here 
that at present a complete evaluation of nutritional quality is im- 
possible without feeding trials with animals. In many instances, 
however, it is not feasible, for ecomonic or other reasons, to subject 
experimentally grown plant material to animal tests. Chemical an- 
alyses in such cases are of value as indicators of nutritional quality. 
Much of the evidence available today indicates that a-lower crude 
fiber and a higher nitrogen content both favorably influence nutri- 
tional quality. 

From the genetic standpoint, the diploid and tetraploid cultures of 
stock A were strictly comparable in that the tetraploids originated by 
direct chromosome doubling of an inbred line of diploid corn that had 
been self-pollinated for many generations and was therefore essentially 
homozygous. The diploid and tetraploid cultures of the other stocks 
were not as strictly comparable as those of stock A. The tetraploid 
lines that produced the parthenogenetic diploids with which these 
other comparisons were made had not been inbred sufficiently to 
render homozygous all of the genes that may have influenced the 
chemical composition of the plants. But within most of these stocks 
more than a single diploid-tetraploid comparison was made, and it 
is believed that a sufficient number of stocks were involved in the 
experiments to furnish significant results, especially as the observed 
differences were consistent. With the exception of one of the seven 
stocks used in the protein analyses of the stover, the results exhibited 
trends which were all in the same direction for both nitrogen and 
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cellulose content. Analyses made in different seasons also exhibited 
similar trends. 

The observed similarities and differences between the diploid and 
tetraploid corn can be explained on a genetic basis simply by assuming 
that certain genes affecting the chemical composition of the corn 
plant act in a cumulative manner and others do not. The genes 
which control nitrogen accumulation in the grain and in the stover 
apparently exhibited a cumulative effect in that the tetraploids, which 
have four doses of these genes in all parts of the plant other than the 
endosperm, and six doses of the endosperm genes, were consistently 
higher in nitrogen content than the related diploids, which have only 
two doses of these genes in the parts of the plant other than the endo- 
sperm and three doses in the endosperm tissue. The genes regulating 
the synthesis of the fatlike constituents of the ether extract and the 
accumulation of inorganic matter in the ash exhibited no dosage 
effect. Both the diploids and the tetraploids contained essentially 
similar amounts of these constituents. 

The negative effect of chromosome doubling which was noted 
with respect to cellulose content of the cornstover was not unexpected. 
Percentage increases in certain constituents must be balanced by 
percentage decreases in other constituents. The decrease in the 
amount of cellulose corresponded roughly with the increase in crude 
protein (N 6.25) content of the tetraploid. 

Randolph and Hand (10) accounted for an increase of 40 percent 
in the content of the carotenoid pigments in yellow corn meal on the 
assumption that the genes for yellow endosperm exerted a cumulative 
effect following chromosome doubling. The present investigation 
has demonstrated a similar but less pronounced increase in the nitrogen 
content of the corn grain and stover, which is interpreted in the same 
way. 

Recently Barr and Newcomer (2) reported that tetraploid cabbage 
contains 36.48 percent more sugars, 23.86 percent more ascorbic acid, 
and 32.62 percent more colloidal nitrogen than diploid cabbage; but 
the diploid had about 14 percent more soluble nitrogen than the 
tetraploid. In experiments with autotetraploid Nicotiana, Noguti, 
Okuma, and Oka (7) and Noguti, Oka, and Otuka (6) reported that 
the tetraploids contained from 51 to 138 percent more nicotine than 
the corresponding diploids and that the content of nitrates, organic 
acids, sugars, total nitrogen, and ether extracts also was increased in 
the tetraploids. 

In certain plants chromosome doubling apparently does not produce 
significant changes in chemical composition. For example, Kostoff 
and Axamitnaja (5) found that diploid and tetraploid petunias did 
not differ appreciably in the amounts of the more important chemical 
constituents. However, tetraploid tomatoes were reported by 
these workers to contain appreciably more water, nitrogen, proteins, 
and ash, but less carbohydrates than the related diploids. But the 
chemical composition of tetraploid tomatoes examined by Fabergé 
(4) was similar to that of the diploids. 

If the chemical composition of plant tissues is regulated by the 
action of specific genes that may or may not act in a cumulative 
manner, it would be expected that the autotetraploids of certain 
species and of certain genotypes within species might differ from those 
of other species or genotypes with respect to chemical composition, 
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and others might not. Differences of this sort were noted in the 
present study and have been reported by other workers. Further 
investigations are needed to determine the general significance of 
chromosome doubling in relation to the nutritional value of culti- 
vated plants. 

SUMMARY AND CONCLUSION 


In certain stocks of corn (Zea maysl.), doubling the number of 
chromosomes was accompanied by an increase in the amount of nitro- 
gen present in the grain and in the stover, the average increase being 
15 percent in the grain and from 20 to 34 percent in the stover. 

The chemical composition of the grain in diploid and tetraploid 
maize was essentially the same with respect to ether extract, ash, 
and crude fiber. 

The diploid and tetraploid corn stover contained essentially the 
same relative amounts of ash, ether extract, and lignin. 

The higher crude protein (N X 6.25) content of the tetraploid 
stover was accompanied by a correspondingly lower content of crude 
fiber and cellulose. 

The changes in chemical composition that accompanied the trans- 
formation from diploidy to tetraploidy are assumed to be due to the 
cumulative action of certain genes, notably those concerned with 
protein metabolism, and the noncumulative action of other genes. 
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THE RELATION OF TEST WEIGHT AND PROTEIN CON- 
TENT TO THE MILLING AND BAKING QUALITY OF 
HARD RED SPRING WHEAT HYBRIDS* 


By R. H. Harris, cereal technologist, and L. R. WALDRON, plant breeder, North 
Dakota Agricultural Experiment Station’ 


INTRODUCTION 


For several years data have accumulated at this station on the micro 
milling and baking quality of a large number of selections from hard 
red spring wheat “hybrids. These “hybrids originated in the wheat 
nursery at Fargo and have been submitted for technologic tests because 
of agronomic promise. The purpose of these tests is to obtain infor- 
mation on their suitability for commercial flour production. While 
many are found to be undesirable, as they show no improvement over 

varieties already grown commercially, a 5 are discovered each year 
which justify further study. If some of these wheats are proved to be 
equal or superior in quality to those already in production they may be 
included in plot trials at one or more localities in the State, and sam- 
ples from these plets are then milled on the Allis-Chalmers experi- 
mental mill and the resultant flours baked by the 100-gm. method. 

The present inquiry was undertaken to discover whether reliable 
preliminary information regarding milling and baking quality, as 
assessed by micro methods, could be secured from a knowledge of test 
weight per bushel and w heat protein content. The situation encoun- 
tered in small-size nursery tests is quite different from the more com- 
mon one of evaluating a number of varied samples consisting of only 
a few wheat varieties, and the relations between test weight and flour 
yield, and wheat protein content and loaf volume, existing in these 
samples might not be found in nursery mater ial. 

Comparisons between the standard and micro milling and baking 
methods reported in the literature increase confidence in micro 
methods. Geddes and Aitken (3), § Geddes and Frisell (4), Harris 
and Sanderson (5), and Sibbitt, Se ott, and Harris (8) hae presented 
the merits of micro procedures and have compared the data with those 
obtained by the use of the Allis-Chalmers mill and 100-gm. baking 
method. 

In the present paper the data secured from four successive wheat 

are discussed, and the information derived from micro milling 
and baking techniques in plant breeding work is pointed out. 


~ Received for publication August 29, 1944. 

*The authors wish to acknowledge the valuable assistance of L. D. Sibbitt in 
performing the milling and baking tests, and of Muriel Elledge in analyzing the 
data. 

* Italic numbers in parentheses refer to Literature Cited, p. 135. 
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EXPERIMENTAL MATERIAL AND METHODS 


The wheats studied were grown in 1940, 1941, 1942, and 1943 in 
the nursery rod rows at Fargo, N. Dak. They were all hard red spring 
wheats of diverse parentage. The individual wheat samples were 
hand-harvested and carefully protected from unfavorable weather 
while in the shock. They were threshed in a small, motor-driven 
thresher when dry and in suitable condition. There were 145 samples 
in 1940 and 1941, 214 in 1942, and 155 in 1943. The samples for tech- 
nologic testing were each composited from the seed from several repli- 
cations. Usually 250 gm. were received for milling and baking tests 
and nitrogen determinations. 

The wheats were cleaned and tempered before milling. The milling 
and baking methods employed were those described by Sibbitt, Scott, 
and Harris (8). Test weight was ascertained by a micro method re- 
ported by Harris and Sibbitt (6) since the quantity of wheat was too 
small for the standard quart procedure, which requires approximately 
114 quarts of grain. Total wheat protein was determined by the 
Kjeldahl procedure. Flour yield is expressed as percent of wheat at 
13.5-percent moisture content. Wheat protein was calculated on a 13.5- 
percent moisture basis, while loaf volumes are from 25 gm. of flour 
containing 86.5 percent dry matter. The doughs were mixed in a Ho- 
bart C-10 mixer equipped with special dough hooks for mixing 50-gm. 
doughs. After mixing, the doughs were divided into equal portions by 
weight and fermented for 3 hours. 


EXPERIMENTAL RESULTS 


The individual determinations were too numerous for presentation, 
and are in themselves of no particular interest. The discussion of the 
results will therefore be limited to the means, the variability, and the 
relations between suitable pairs of the four variables: test weight per 
bushel, flour yield, wheat protein content, and loaf volume. In this 
discussion differences significant at the 5-percent point will be defined 
as significant, while those significant at the 1-percent point will be 
called very significant. 


MEANS 


The mean values for each variable studied for each year are shown 
in table 1. There are very significant differences among these annual 
values. The differences in test weight between 1940 and 1948 and be- 
tween 1941 and 1942 are very significant, while those between 1940 
and 1942 and between 1941 and 1943, are significant at the 5-percent 
point. These differences are principally caused by weather condi- 
tions, no doubt, although the experimental material varies in genetic 
composition from year to year. For flour yield all yearly averages 
differ very significantly, except those for 1942 and 1948, which are 
identical. There seems to be little association bet ween the correspond- 
ing means for test weight and flour yield, despite a priori expectations. 
For wheat protein very significant differences exist between 1940 and 
1941 and between 1941 and 1942, while the 1943 values are very 
significantly higher than any of the other three. The difference be- 
tween the 1940 and 1942 means is just significant. There is possibly 
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some degree of association between the yearly averages for test weight 
and wheat protein since the 2 years of high test weight were those in 
which the protein was lowest, while low test weight crops apparently 
had high protein content. All loaf volume differences were very 
significant. Little relation between the wheat protein and loaf vol- 
ume means are evident. 


STANDARD DEVIATIONS AND COEFFICIENTS OF VARIATION 


The standard deviations and coefficients of variation for each vari- 
able are shown in table 1 for each year. The significance of differences 
in variability between years was not tested, but there seems to be no 
extreme differences. Apparently variations in weather conditions and 
the genetic composition of the wheats did not greatly affect the vari- 
ability of the results. 


TABLE 1.—Means, standard deviations, and coefficients of variation of test weight, 
flour yield, wheat protein, and loaf volume of wheats tested for each year. 
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MAXIMUM AND MINIMUM VALUES 


The maximum and minimum values of wheats tested, with corre- 
sponding ranges, are shown in table 2. It is realized that these val- 
ues do not reflect the variability of the material as accurately as stand- 
ard deviations since they are determined by very few variates and 
are subject to large errors of random sampling. They are presented 
merely to show the substantial differences in actual values found in 
a single crop year, and to illustrate the fact that absence of signifi- 
‘ant correlation between suitable pairs of variables was not caused 
by lack of variation in the create 5 which was all grown in the same 
location. Test weight varied 9.4 pounds per bushel in 1943, and flour 
yield varied 12.6 percent for the same crop. In 1942 wheat protein 
content differed 5.8 percent and loaf volume 120 ce. 

The causes of these differences within years probably lie in varia- 
tions in the genetic composition of the wheats, causing them to respond 
differently to weather conditions, since the soil properties should be 
strictly comparable for all samples. 


CORRELATION COEFFICIENTS 


The correlation coefficients between test weight and flour yield, and 
between wheat protein and loaf volume for each year are shown in 
table 3. Little comment on these is required, except to point out the 
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rather surprising lack of close association in most instances between 
the different pairs of variables. With the exception of wheat protein 
und loaf volumes for 1941 and 1948 these are too small to be of much 
practical significance. For instance, no worth while information 
regarding flour yield can be secured from a knowledge of test weight, 
and the only information of value from the test weight data is for 
grading and general descriptive purposes. The correlations between 
wheat protein and loaf volume are considerably better, but for the 
1940 crop there is no relation of statistical significance, while for the 
1942 crop the correlation, though significant, is of little practical 
ralue. For the remaining 2 years, 1941 and 1943, the correlations are 
relatively high but not high enough to permit the prediction of loaf 
volume from wheat protein with any reasonable degree of confidence. 
The relations between some of these associated data were examined 
for nonlinearity but no evidence of this was found. 


TABLE 2.—Mazximum and minimum values of wheats tested with corresponding 


ranges of test weight, wheat protein, flour yield, and loaf volume for each 
year. 
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TABLE 3.—Correlation cofficients between test weight and flour yield, and between 
wheat protein and loaf volume for each year 


Variables correlated | Crop Correlation 





i | year = 
_—-— -- ——. | — ee SesS ae enna retes 
| | 1940 +0. 108 
Test weight (pounds per bushel) _.__-.-- | Flour yield (percent) _..............-...-.. — = 4 
| \{ 1943 +. 136 
| | 1940 +. 076 
Wheat protein (percent) -.........---..- Loaf volume (cubic centimeter) __._______- | on > 4 
| | 1943 +. 536 





1 The value of r,, at-the 5-percent point for the different years was as follows: 1940 and 1941, 0.163; 1942, 


0.134; 1943, 0.158. Its value at the 1-percent point was: 1940 and 1941, 0.215; 1942, 0.177; 1943, 0.208. 
DISCUSSION 


The absence of high correlation between the two principal pairs of 
variables shows very clearly the fallacy of any but the most general 
and uncertain prediction of flour yield and loaf volume from test 
weight and wheat protein content, respectively. This is particularly 
true respecting test weight and flour yield. The correlation coeffi- 
cients between wheat protein content and loaf volume for the 4 years 
were transformed and combined by the method suggested by Fisher 
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(2). The value obtained for the combined results was only +0.370. 
Since the square of this value is 0.137, one could conclude that wheat 
protein content accounted for only approximately 14 percent of the 
loaf volume variation, indicating that about 86 percent is dependent 
upon other factors. The variation in flour yield independent of test 
weight would be greater because the seasonal correlations between 
these two variables are lower than for wheat protein content and loaf 
volume. Thus the plant breeder basing his estimates of the flour- 
yielding and baking capacity of his hybrids and selections upon the 
relatively easy and rapid determinations of test weight and protein 
content would be subject to substantial errors in his decisions regard- 
ing the milling and baking quality of his wheats, as determined by 
micro milling and baking methods, with the expenditure of needless 
time and expense in later eradicating inferior selections from the 
breeding program. In the interim there would be excellent prospects 
of some really promising varieties being discarded. 

The results obtained from this investigation conducted with micro 
methods are in essential agreement with those secured by Hayes et al. 
(7) and Ausemus et al. (7) by the use of the Allis-Chalmers experi- 
mental mill and the 100-gm. baking method on wheats of diverse 
genetic origin grown at four Minnesota stations. However, Waldron 
and Mangels (9), also employing the Allis-Chalmers mill and the 
100-gm. baking method, found significant relationships between flour 
protein content and loaf volume for nursery tests over a 4-year period 
(r= +0441). Their values were not high enough, though, to be 
useful for predicting loaf volume from flour protein content. 


SUMMARY AND CONCLUSIONS 


Very significant differences between crop years were found in test 
weight, wheat protein content, flour yield, and flour loaf volume for 
659 samples of hard red spring wheat produced in the Fargo nursery 
over a 4-year period. These differences were no doubt caused princi- 
pally by yearly variation in weather conditions, with variations in 
the genetic composition of the wheats being less important. 

Correlation coefficients for each of the four crop years between test 
weight and flour yield were so low as to be of no practical significance, 
while the coefficients between wheat protein content and loaf volume, 
though generally higher, were of little practical value for predicting 
loaf volume from wheat protein content. The relationships between 
the different pairs of variables were essentially linear in character. 

The results secured in these investigations show the danger of rely- 
ing solely upon test weight and wheat protein content data to predict 
accurately milling and baking quality in nursery rod-row samples of 
hard red spring wheat. 
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